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Abstract I give a unified perspective on the properties of a variety of quantum liq- 
uids using the theory of quantum phase transitions. A central role is played by a zero 
density quantum critical point which is argued to control the properties of the dilute 
gas. An exact renormalization group analysis of such quantum critical points leads 
to a computation of the universal properties of the dilute Bose gas and the spinful 
Fermi gas near a Feshbach resonance. 



1 Introduction 

This article is adapted from Chapter 16 of \Quantum Phase Transitions, 2nd edition| 
Cambridge University Press. 

It is not conventional to think of dilute quantum liquids as being in the vicin- 
ity of a quantum phase transition. However, there is a simple sense in which they 
are, although there is often no broken symmetry or order parameter associated with 
this quantum phase transition. We shall show below that the perspective of such a 
quantum phase transition allows a unified and efficient description of the universal 
properties of quantum liquids. 

Stated most generally, consider a quantum liquid with a global U(l) symmetry. 
We shall be particularly interested in the behavior of the conserved density, gener- 
ically denoted by Q (usually the particle number), associated with this symmetry. 
The quantum phase transition is between two phases with a specific T = behav- 
ior in the expectation value of Q. In one of the phases, {Q) is pinned precisely at a 
quantized value (often zero) and does not vary as microscopic parameters are var- 
ied. This quantization ends at the quantum critical point with a discontinuity in the 
derivative of {Q) with respect to the tuning parameter (usually the chemical poten- 
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tial), and {Q) varies smoothly in the other phase; there is no discontinuity in the 
value of {Q), however. 

The most familiar model exhibiting such a quantum phase transition is the dilute 
Bose gas. We express its coherent state partition function, Zb, in terms of complex 
field ^^(x, t), where jc is a t/-dimensional spatial co-ordinate and T is imaginary 
time: 

Zb = y WB(x,T)exp j^J^ dx J d''x^B 

J^s = f; ^+ ^IVftp-Ml^l' + ylffer (1) 
We can identify the charge Q with the boson density "fg "ffe 

{Q) = -^ = m'), (2) 

with = ^ (r /y ) InZfi. The quantum critical point is precisely at /i = and T = 0, 
and there are no fluctuation corrections to this location from the terms in So at 
T ^0, {Q) takes the quantized value (g) = for jU < 0, and (g) > for jU > 0; we 
will describe the nature of the onset at = and finite-T crossovers in its vicinity. 

Actually, we will begin our analysis in Section [2] by a model simpler than Zb, 
which displays a quantum phase transition with the same behavior in a conserved 
U(l) density (Q) and has many similarities in its physical properties. The model is 
exactly solvable and is expressed in terms of a continuum canonical spinless fermion 
field f^; its partition function is 

Zf = y WF(x,T)exp ^- J^^^ dr J d'^x^F 

J^F = Wp ^+ ±|Vff|2-^|ff|2. (3) 
oT Zm 

J^ff is just a free field theory. Like Zg, Zf has a quantum critical point at —0, 
T = and we will discuss its properties; in particular, we will show that all possible 
fermionic nonlinearities are irrelevant near it. The reader should not be misled by 
the apparently trivial nature of the model in (|3); using the theory of quantum phase 
transitions to understand free fermions might seem like technological overkill. We 
will see that Zp exhibits crossovers that are quite similar to those near far more 
complicated quantum critical points, and observing them in this simple context leads 
to considerable insight. 

In general spatial dimension, d, the continuum theories Zb and Zf have differ- 
ent, though closely related, universal properties. However, we will argue that the 
quantum critical points of these theories are exactly equivalent ind = I. We will see 
that the bosonic theory Zb is strongly coupled in t/ = 1, and will note compelling 
evidence that the solvable fermionic theory Zf is its exactly universal solution in 
the vicinity of the /i = 0, T = quantum critical point. This equivalence extends to 
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observable operators in both theories, and allows exact computation of a number of 
universal properties of Zgind ^ I. 

Our last main topic will be a discussion of the dilute spinful Fermi gas in Sec- 
tion|4] This generalizes Zp to a spin 5=1/2 fermion ^pa, with <T =t, i- Now Fermi 
statistics do allow a contact quartic interaction, and so we have 

Zf , = J &Wf^ {x, x) &Wfi {x,T)exp(^- J^^'^ dr J d''x 

= f;^ ^ + ^ I vff al' - Ml-iVal' + Mofn'f'n^'i'n- (4) 

This theory conserves fermion number, and has a phase transition as a function of 
increasing /i from a state with fermion number to a state with non-zero fermion 
density. However, unlike the above two cases of Zb and Zf, the transition is not 
always at = 0. The problem defined in (|4|i has recently found remarkable ex- 
perimental applications in the study of ultracold gases of fermionic atoms. These 
experiments are aslo able to tune the value of the interaction mq over a wide range of 
values, extended from repulsive to attractive. For the attractive case, the two-particle 
scattering amplitude has a Feshbach resonance where the scattering length diverges, 
and we obtain the unitarity limit. We will see that this Feshbach resonance plays a 
crucial role in the phase transition obtained by changing /i, and leads to a rich phase 
diagram of the so-called "unitary Fermi gas". 

Our treatment of Zfs in the experimental important case of d ~ 3> will show 
that it defines a strongly coupled field theory in the vicinity of the Feshbach reso- 
nance for attractive interactions. It therefore pays to find alternative formulations of 
this regime of the unitary Fermi gas. One powerful approach is to promote the two 
fermion bound state to a separate canonical Bose field. This yields a model, Zfb 
with both elementary fermions and bosons ; i.e. it is a combination of Zg and Zfs 
with interactions between the fermions and bosons. We will define Zfb in SectionH) 
and use it to obtain a number of experimentally relevant results for the unitary Fermi 
gas. 

Section 12] will present a thorough discussion of the universal properties of Z/r. 
This will be followed by an analysis of Zb in Section [3] where we will use renor- 
malization group methods to obtain perturbative predictions for universal properties. 
The spinful Fermi gas will be discussed in Section |4l 



2 The Dilute Spinless Fermi Gas 

This section will study the properties of Zf in the vicinity of its jU = 0, 7 = 
quantum critical point. As Zf is a simple free field theory, all results can be obtained 
exactly and are not particularly profound in themselves. Our main purpose is to show 
how the results are interpreted in a scaling perspective and to obtain general lessons 
on the nature of crossovers at T > 0. 
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First, let us review the basic nature of the quantum critical point at T = 0. A 
useful diagnostic for this is the conserved density Q, which in the present model 
we identify as f^f^. As a function of the tuning parameter /i, this quantity has a 
critical singularity at ;U = 0: 

(u,^u,\^HSdld){2miiYI^, ju>0, 
\^F^P/~\ 0, M<0, 

where the phase space factor Sd = 2/[r{d /2)(4;r)'^'/^]. 

We now proceed to a scaling analysis. Notice that at the quantum critical point 
/i = 0, r = 0, the theory is invariant under the scaling transformations: 

/ -I 
X = xe , 

t' = Tf-^'^ (6) 

provided we make the choice of the dynamic exponent 

z = 2. (7) 

The parameter m is assumed to remain invariant under the rescaling, and its role is 
simply to ensure that the relative physical dimensions of space and time are com- 
patible. The transformation (|6) also identifies the scaling dimension 

dim[fF] = d/2. (8) 

Now turning on a nonzero jj., it is easy to see that is a relevant perturbation 
with 

dim[^] = 2. (9) 

There will be no other relevant perturbations at this quantum critical point, and so 
we have for the correlation length exponent 

v = l/2. (10) 

We can now examine the consequences of adding interactions to ^p. A contact 
interaction such as / dx{^p{x)^f{x))^ vanishes because of the fermion anticommu- 
tation relation. (A contact interaction is however permitted for a spin- 1/2 Fermi gas 
and will be discussed in Section|4]i The simplest allowed term for the spinless Fermi 
gas is 

i^l =A('P;(.^,T)V'P;(x,T)fF(.^,T)V'f^F(x,T)), (11) 

where A is a coupling constant measuring the strength of the interaction. However, 
a simple analysis shows that 

dim[A] = -d. (12) 
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This is negative and so X is in^elevant and can be neglected in the computation of 
universal crossovers near the point /i = T = 0. In particular, it will modify the result 
(|5]l only by contributions that are higher order in /i. 

Turning to nonzero temperatures, we can write down scaling forms. Let us define 
the fermion Green's function 

GF{x,t) = {WF{x,t)W^p{Q,Q))- (13) 

then the scaUng dimensions above imply that it satisfies 

Gf (x,f) = {2mTYI^<pG, ((2mr)i/2x,rf,^) , (14) 

where is a fully universal scaling function. For this particularly simple theory 
we can of course obtain the result for Gf in closed form: 

f d'^k e'^-'(*V(2'«)-M)' 

and it is easy to verify that this obeys the scaling form (fl4] i. Similarly the free energy 
^F has scaling dimension c/ + z, and we have 

^, = T"I^+'^,^J^\ (16) 



with a. universal scaling function; the explicit result is, of course, 

(2 

which clearly obeys ( fTSI l. The crossover behavior of the fermion density 



{Q) = {%^f)=—^ (18) 

follows by taking the appropriate derivative of the free energy. Examination of these 
results leads to the crossover phase diagram of Fig. [T] We will examine each of 
the regions of the phase diagram in turn, beginning with the two low-temperature 
regions. 



2.1 Dilute Classical Gas, /c^r ^ /i < 

The ground state for /i < is the vacuum with no particles. Turning on a nonzero 
temperature produces particles with a small nonzero density ^e^l''!/^. The de 
Broglie wavelength of the particles is of order T^'/^, which is significantly smaller 
than the mean spacing between the particles, which diverges as g!''!/'^'^ as T 0. 
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Fig. 1 Phase diagram of the dilute Fermi gas Zf (Eqn. {3j) as a function of the chemical potential 
^ and the temperature T . The regions are separated by crossovers denoted by dashed lines, and 
their physical properties are discussed in the text. The full lines are contours of equal density, 
with higher densities above lower densities; the zero density line is < 0, T = 0. The line ;U > 0, 
r = is a line of z = 1 critical points that controls the longest scale properties of the low-T Fermi 
liquid region. The critical end point /i = 0, T = has z = 2 and controls global structure of the 
phase diagram. In d = 1, the Fernii liquid is more appropriately labeled a Tomonaga-Luttinger 
liquid. The shaded region marks the boundary of applicability of the continuum theory and occurs 
at ii,T r^w. 

This implies that the particles behave semiclassically. To leading order from (IT5i . 
the fermion Green's function is simply the Feynman propagator of a single particle 



and the exclusion of states from the other particles has only an exponentially small 
effect. Notice that Gp is independent of n and T and (fT9] l is the exact result for 
jj. = T = 0. The free energy, from (fT&l and dTTi . is that of a classical Boltzmann gas 





(20) 



2.2 Fermi Liquid, kgT ^ /i, /i > 



The behavior in this regime is quite complex and rich. As we will see, and as noted in 
Fig.[T] the line jj. > 0, T — is itself a line of quantum critical points. The interplay 
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between these critical points and those of the ^ = 0, T = critical end point is 
displayed quite instructively in the exact results for Gf and is worth examining in 
detail. It must be noted that the scaling dimensions and critical exponents of these 
two sets of critical points need not (and indeed will not) be the same. The behavior 
of the > 0, r = critical line emerges as a particular scaling limit of the global 
scaling functions of the = 0, T = critical end point. Thus the latter scaling 
functions are globally valid everywhere in Fig.[T] and describe the physics of all its 
regimes. 

First it can be argued, for example, by studying asymptotics of the integral in 
(fTSl ), that for very short times or distances, the correlators do not notice the conse- 
quences of other particles present because of a nonzero T or ji and are therefore 
given by the single-particle propagator, which is the T = = result in (fT9] l. More 
precisely we have 

G(Jc,f)isgivenby([T9llfor|x|<(2m^)"'/^ kl < -• (21) 

With increasing x or f, the restrictions in (ISTT i are eventually violated and the con- 
sequences of the presence of other particles, resulting from a nonzero jx, become 
apparent. Notice that because jj. is much larger than T, it is the first energy scale to 
be noticed, and as a first approximation to understand the behavior at larger x we 
may ignore the effects of T. 

Let us therefore discuss the ground state for /i > 0. It consists of a filled Fermi sea 
of particles (a Fermi liquid) with momenta k < kf ^ (Imjiy^^. An important prop- 
erty of the this state is that it permits excitations at arbitrarily low energies (i.e., it is 
gapless). These low energy excitations correspond to changes in occupation number 
of fermions arbitrarily close to As a consequence of these gapless excitations, 
the points > (T = 0) form a line of quantum critical points, as claimed earlier. 
We will now derive the continuum field theory associated with this line of critical 
points. We are interested here only in x and t values that violate the constraints in 
(l2Tl l. and so in occupation of states with momenta near ±kF- So let us parameterize, 
in fif = 1, 

^{x, t) ^ e'^^'^Rix, t) + e-^'^-^^fi (jc, t), (22) 

where Wr l describe right- and left-moving fermions and are fields that vary slowly 
on spatial scales ^l/kp = (l/2m/i)'/^ and temporal scales most of the 

results discussed below hold, with small modifications, in all d. Inserting the above 
parameterization in and keeping only terms lowest order in spatial gradients, 
we obtain the "effective" Lagrangean for the Fermi liquid region, J/fFi in d = 1 : 

where vp ^ kF / m = [Ipi / niY 1^ is the Fermi velocity. Now notice that ^fl is invari- 
ant under a scaling transformation, which is rather different from (|6]l for the jJ. = 0, 
T = quantum critical point: 
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X' = 



xe 



t 



%^^{x',x')^WR,i^{x,x)e'l\ 



(24) 




(25) 



Aim[WR,L] = 1/2, 



(26) 



which actually holds for all d and therefore differs from (|8). Further notice that Vf, 
and therefore /i, are invariant under rescaling, unlike (|9) at the /i = critical point. 
Thus VF plays a role rather analogous to that of m at the ji = critical point: It 
is simply the physical units of spatial and length scales. The transformations (l24t 
show that is scale invariant for each value of ji, and we therefore have a line 
of quantum critical points as claimed earlier. It should also be emphasized that the 
scaling dimension of interactions such as A will also change; in particular not all 
interactions are irrelevant about the fi ^ critical points. These new interactions 
are, however, small in magnitude provided jU is small (i.e., provided we are within 
the domain of validity of the global scaling forms (fl4] i and (fTSI l, and so we will 
neglect them here. Their main consequence is to change the scaling dimension of 
certain operators, but they preserve the relativistic and conformal invariance of ^fl- 
This more general theory of d ~ I fermions is the Tomonaga-Luttinger liquid. 



2.3 High-T Limit, kgl 

This is the last, and in many ways the most interesting, region of Fig.[T] Now T is 
the most important energy scale controlling the deviation from the jj. = 0, T = 
quantum critical point, and the properties will therefore have some similarities to 
the "quantum critical region" of other strongly interacting models ifTTl . It should be 
emphasized that while the value of T is significantly larger than |jLl |, it cannot be so 
large that it exceeds the limits of applicability for the continuum action ^f- If we 
imagine that ^f was obtained from a model of lattice fermions with bandwidth w, 
then we must have T -^w. 

We discuss first the behavior of the fermion density. In the high-T limit of the 
continuum theory ^f, |/i | <C T ^ w, we have from (fTTI i and ( fTSl l the universal result 




Dilute Fenni and Bose Gases 



9 



= (2,„rrV2^(^/2)(L_p. (27) 

This density implies an interparticle spacing that is of order the de Broglie wave- 
length = (l/2mr)'/^. Hence thermal and quantum effects are to be equally impor- 
tant, and neither dominate. 

For completeness, let us also consider the fermion density for T ':$>w (the region 
above the shaded region in Fig. [T]i, to illustrate the limitations on the continuum 
description discussed above. Now the result depends upon the details of the nonuni- 
versal fermion dispersion; on a hypercubic lattice with dispersion — /i, we obtain 



'^A' d''k 1 



,1a {2nY e(^k-ii)lT + 1 



1 1 f^'l" d''k 



2a'' AT J-n/a [2% 



{et^li)^&{\/T^). (28) 



The limits on the integration, which extend from — tt/a to n/a for each momentum 
component, had previously been sent to infinity in the continuum limit a ^ 0. In the 
presence of lattice cutoff, we are able to make a naive expansion of the integrand 
in powers of \/T, and the result therefore only contains negative integer powers of 
T . Contrast this with the universal continuum result (|27| where we had noninteger 
powers of T dependent upon the scaling dimension of 'f'. 

We return to the universal high- T region, | | <C T <C w, and describe the behavior 
of the fermionic Green's function Gf, given in ( fTsT l. At the shortest scales we again 
have the free quantum particle behavior of the ^ — Q,T — Q critical point: 

G/r(x,f) is given by ^ for \x\ < {2mTy^l^%\t\ < ^. (29) 

Notice that the limits on x and t in ( |29] l are different from those in (1211 . in that they 
are determined by T and not /i. At larger \x\ or f the presence of the other thermally 
excited particles becomes apparent, and Gf crosses over to a novel behavior charac- 
teristic of the high-T region. We illustrate this by looking at the large-x asymptotics 
of the equal-time G in t/ = 1 (other d are quite similar): 

r dk e'^'^ 

^-(^'«)=y2^TT;^- ^''^ 

For large x this can be evaluated by a contour integration, which picks up contribu- 
tions from the poles at which the denominator vanishes in the complex k plane. The 
dominant contributions come from the poles closest to the real axis, and give the 
leading result 

2 \ 1/2 



Gf{\x\ ^ oo,0) = -[ —,] exp(-(l-/)(m;rr)'/^x). (31) 
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Thermal effects therefore lead to an exponential decay of equal-time correlations, 
with a correlation length t, = (mTzT) ' . Notice that the T dependence is precisely 
that expected from the exponent z = 2 associated with the /i = quantum critical 
point and the general scaling relation ^ ^ T^^I-\ The additional oscillatory term 
in dsTl ) is a reminder that quantum effects are still present at the scale , which is 
clearly of order the de Broglie wavelength of the particles. 



3 The Dilute Bose Gas 

This section will study the universal properties quantum phase transition of the di- 
lute Bose gas model Zb in ([T]) in general dimensions. We will begin with a simple 
scaling analysis that will show that (i = 2 is the upper-critical dimension. The first 
subsection will analyze the case d <2'm some more detail, while the next subsec- 
tion will consider the somewhat different properties in c/ = 3. Some of the results of 
this section were also obtained by Kolomeisky and Straley 191 [TOl . 

We begin with the analog of the simple scaling considerations presented at the 
beginning of Section |2] At the coupling u = 0, the ^ = quantum critical point of 
is invariant under the transformations (|6), after the replacement f^, and 

we have as before z = 2 and 

dim[»f'B] =£//2, dim[^]=2; (32) 

these results will shortly be seen to be exact in all d. We can easily determine the 
scaling dimension of the quartic coupling u at the m = 0, /i = fixed point under the 
bosonic analog of the transformations (|6); we find 

dim[Mo] = 2~d. (33) 

Thus the free-field fixed point is stable for li > 2, in which case it is suspected that 
a simple perturbative analysis of the consequences of u will be adequate. However, 
for d <2,?L more careful renormalization group-based resummation of the conse- 
quences of u is required. This identifies (i = 2 as the upper-critical dimension of the 
present quantum critical point. 

Our analysis of the case d <2 for the dilute Bose gas quantum critical point will 
find, somewhat surprisingly, that all the renormalizations, and the associated flow 
equations, can be determined exactly in closed form. We begin by considering the 
one-loop renormalization of the quartic coupling mq at the /i = 0, T = quantum 
critical point. It turns out that only the ladder series of Feynman diagrams shown in 
Fig. |2] need be considered (the T matrix). Evaluating the first term of the series in 
Fig.|2]for the case of zero external frequency and momenta, we obtain the contribu- 
tion 

dco f d''k 1 1 , /■ d^'k m 



dco r 

27t J 



{2nY {-ico + k^/{2m)) {ico + k^/{2m)) V {2ny k^ 
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-1 — * I * -"^n — * I * — r 

I I I I I 



n — ' I — — r 



+ ' ' ' ' + ... 

Fig. 2 The ladder series of diagrams that contribute the renormalization of the coupling u in Zb 
for d<2. 



(the remaining ladder diagrams are powers of (|34) and form a simple geometric 
series). Notice the infrared singularity for d <2, which is cured by moving away 
from the quantum critical point, or by external momenta. 

We can proceed further by a simple application of the momentum shell RG. Note 
that we will apply cutoff A only in momentum space. The RG then proceeds by 
integrating all frequencies, and momentum modes in the shell between Ae^' and 
A. The renormalization of the coupling uq is then given by the first diagram in 
Fig.|2] and after absorbing some phase space factors by a redefinition of interaction 
coupling 

MO = i^—^u, (35) 

we obtain ||6l U\ 

du 

dl-'^-Y- ^''^ 

Here ~ 2/{r{d /2)(47r)'^/^) is the usual phase space factor, and 

e = 2-d. (37) 

Note that for e > 0, there is a stable fixed point at 

M*==2£, (38) 

which will control all the universal properties of Zb- 

The flow equation (|36] |, and the fixed point value (l38T l are exact to all orders 
in II or e, and it is not necessary to consider iz-dependent renormalizations to the 
field scale of "ffe or any of the other couplings in Zb- This result is ultimately a 
consequence of a very simple fact: The ground state of Zg at the quantum critical 
point jLi = is simply the empty vacuum with no particles. So any interactions that 
appear are entirely due to particles that have been created by the external fields. In 
particular, if we introduce the bosonic Green's function (the analog of (flSl l) 

GB(x,0 = (fB(x,f)f;(0,0)), (39) 
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then for jJ. <Q and T — 0, its Fourier transform G{k, (o) is given exactly by the free 
field expression 

GB{k,co) = , . . (40) 

The field creates a particle that travels freely until its annihilation at (x,f) by the 
field "ffe; there are no other particles present at T = 0, /i < 0, and so the propagator 
is just the free field one. The simple result (|40] | implies that the scaling dimensions 
in ( l32l ) are exact. Turning to the renormalization of «, it is clear from the diagram 
in Fig. ID that we are considering the interactions of just two particles. For these, 
the only nonzero diagrams are the one shown in Fig. |2] which involve repeated 
scattering of just these particles. Formally, it is possible to write down many other 
diagrams that could contribute to the renormalization of u; however, all of these 
vanish upon performing the integral over internal frequencies for there is always one 
integral that can be closed in one half of the frequency plane where the integrand 
has no poles. This absence of poles is of course just a more mathematical way of 
stating that there are no other particles around. 

We will consider application of these renormalization group results separately 
for the cases below and above the upper-critical dimension of c/ = 2. 



3.1 d<2 

First, let us note some important general implications of the theory controlled by 
the fixed point interaction ( |38] |. As we have already noted, the scaling dimensions 
of and n are given precisely by their free field values in d32t , and the dynamic 
exponent z also retains the tree-level value z = 2. All these scaling dimensions are 
identical to those obtained for the case of the spinless Fermi gas in Section |2] Fur- 
ther, the presence of a nonzero and universal interaction strength u* in ( |38] | implies 
that the bosonic system is stable for the case ji > because the repulsive interac- 
tions will prevent the condensation of infinite density of bosons (no such interaction 
was necessary for the fermion case, as the Pauli exclusion was already sufficient to 
stabilize the system). These two facts imply that the formal scaling structure of the 
bosonic fixed point being considered here is identical to that of the fermionic one 
considered in Section|2]and that the scaling forms of the two theories are identical. 
In particular, Gb will obey a scaling form identical to that for Of in (fl4] | (with a cor- 
responding scaling function ^Gb)^ while the free energy, and associated derivatives, 
obey ( fT6b (with a scaling function 'P,^g)- The universal functions <l>Qg and <Pi^g can 
be determined order by order in the present e =2 — d expansion, and this will be 
illustrated shortly. 

Although the fermionic and bosonic fixed points share the same scaling dimen- 
sions, they are distinct fixed points for general d < 2. However, these two fixed 
points are identical precisely in = 1 ||T9l . Evidence for this was presented in 
Ref. in, where the anomalous dimension of the composite operator was com- 
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puted exactly in the e expansion and was found to be identical to that of the cor- 
responding fermionic operator Assuming the identity of the fixed points, we can 
then make a stronger statement about the universal scaling function: those for the 
free energy (and all its derivatives) are identical = 0^^, ind ~ 1 . In particular, 
from (Tl\ and ( fTSl l we conclude that the boson density is given by 

(g)-W-ft>-/g ^,,„,„,|!„„,^, (41) 

in c/ = 1 only. The operators "ffe and are still distinct and so there is no reason for 
the scaling functions of their correlators to be the same. However, in t/ = 1, we can 
relate the universal scaling function of to those of via a continuum version of 
the Jordan- Wigner transformation 

^B{x,t) = exp (in f dyWXy,t)WF{y,t)] ff (^,f). (42) 



This identity is applied to obtain numerous exact results in Ref. fTTll 

As not all observables can be computed exactly in c/ = 1 by the mapping to the 
free fermions, we will now consider the e = 2 — d expansion. We will present a 
simple e expansion calculation ifTSl for illustrative purposes. We focus on density 
of bosons at r = 0. Knowing that the free energy obeys the analog of ( fTSl l. we can 
conclude that a relationship like Q holds: 

at r = 0, with %i a universal number The identity of the bosonic and fermionic 
theories ind =1 implies from Q or from (l4ll that '^i=Si/l = l/n. We will show 
how to compute %i in the e expansion; similar techniques can be used for almost 
any observable. 

Even though the position of the fixed point is known exactly in ( |38] ). not all 
observables can be computed exactly because they have contributions to arbitrary 
order in u. However, universal results can be obtained order-by-order in m, which 
then become a power series in £ = 2 — (i. As an example, let us examine the low 
order contributions to the boson density. To compute the boson density for jJ. > 0, 
we anticipate that there is condensate of the boson field "f^, and so we write 

VBix,x) = % + '¥iix,t), (44) 

where *Fi has no zero wavevector and frequency component. Inserting this into 
in ([T]), and expanding to second order in fi, we get 

-At |f 1 12 + ^ (4| fo|2 1 fi |2 + ^/jZfr^ + fo*2^/2) . (45) 
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This is a simple quadratic theory in the canonical Bose field 'f'l , and its spectrum and 
ground state energy can be determined by the familiar Bogoliubov transformation. 
Carrying out this step, we obtain the following formal expression for the free energy 
density ^ as a function of the condensate 'fb at T = 0: 




— -M + 2Mo|fo|- 




2m 



-Ai + 2Mo|fo|2 



(46) 



To obtain the physical free energy density, we have to minimize with respect 
to variations in % and to substitute the result back into (|46T l. Finally, we can take 
the derivative of the resulting expression with respect to /i and obtain the required 
expression for the boson density, correct to the first two orders in uq: 



d"k 



y/k^{k^+4-mll) 



(47) 



To convert ( |47] i into a universal result, we need to evaluate it at the coupling 
appropriate to the fixed point (|38] l. This is most easily done by the field-theoretic 
RG. So let us translate the RG equation ( |36] | into this language. We introduce a 
momentum scale p. (the tilde is to prevent confusion with the chemical potential) 
and express uq in terms of a dimensionless coupling by 



(2m)^f UR 



(48) 



The motivation behind the choice of the renormalization factor in (I48b is that the 
renormalized four-point coupling, when expressed in terms of ur, and evaluated 
in = 2 — e, is free of poles in e as can easily be explicitly checked using ( |34] | 
and the associated geometric series. Then, we evaluate (|47t at the fixed point value 
of ur, compute any physical observable as a formal diagrammatic expansion in uq, 
substitute uq in favor of ur using ( |48] l, and expand the resulting expression in powers 
of e. All poles in e should cancel, but the resulting expression will depend upon the 
arbitrary momentum scale /I. At the fixed point value m^, dependence upon /i then 
disappears and a universal answer remains. In this manner we obtain from (|47] | a 
universal expression in the form (|43) with 



1 

2e 



ln2- 1 



(49) 
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Fig. 3 Crossovers of the dilute Bose gas in = 3 as a function of the chemical potential /i and the 
temperature T . The regimes labeled A, B, C are described in Ref. 1171 . The solid line is the finite- 
temperature phase transition where the superfluid order disappears; the shaded region is where 
there is an effective classical description of thermal fluctuations. The contours of constant density 
are similar to those in Fig. [Tjand are not displayed. 



3.2 d = 3 

Now we briefly discuss 2 < d < 4: details appear elsewhere [ITTl . In d = 2, the 
upper critical dimension, there are logarithmic corrections which were computed by 
Prokof 'ev et al. ifTsl . Related results, obtained through somewhat different methods, 
are available in the Hterature llT3l [14] |6] [T9l . 

The quantum critical point at /i = 0, T = is above its upper-critical dimension, 
and we expect mean-field theory to apply. The analog of the mean-field result in the 
present context is the T = relation for the density 

, t \ f M/mqH , M > 0, 

(^;^.)^|^/; ;;<o; '''' 

where the ellipses represents terms that vanish faster as /i — > 0. Notice that this ex- 
pression for the density is not universally dependent upon pi; rather it depends upon 
the strength of the two-body interaction uq (more precisely, it can be related to the 
s-wave scattering length a by uq = Ana/m). The crossovers and phase transitions at 
r > are sketched in Fig. |3] These are similar to those of the spinless Fermi gas, 
but now there can be a phase transition within one of the regions. Explicit expres- 
sions for the crossovers ifTTI have been presented by Rasolt et al. lfT6l . Weichman et 
al. l|29l and also addressed in earher work Ii22. ,23. ,3J . 
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4 The Dilute Spinful Fermi Gas: the Feshbach Resonance 

This section turns to the case of the spinful Fermi gas with short-range interactions; 
as we noted in the introduction, this is a problem which has acquired renewed im- 
portance because of the new experiments on ultracold fermionic atoms. 

The partition function of the theory examined in this section was displayed in 
(lUi. The renormalization group properties of this theory in the zero density limit are 
identical to those the dilute Bose gas considered in Section [3] The scaling dimen- 
sions of the couplings are the same, the scaling dimension of f^tj is d /2 as for ffe 
in ( [32l i, and the flow of the u is given by (|36] l. Thus for d <2, a spinful Fermi gas 
with repulsive interactions is described by the stable fixed point in ( [38] l. 

However, for the case of spinful Fermi gas case, we can consider another regime 
of parameters which is of great experimental importance. We can also allow u to 
be attractive: unlike the Bose gas case, the m < case is not immediately unstable, 
because the Pauli exclusion principle can stabilize a Fermi gas even with attractive 
interactions. Furthermore, at the same time we should also consider the physically 
important case with d >2, when e < 0. The distinct nature of the RG flows predicted 
by ( |36] | for the two signs of e are shown in Fig.|4] 



d<2 





-4- 



U 



(a) 



d>2 





-4- 



(b) 



u 



Fig. 4 The exact RG flow of )36t . (a) For d <2 (e > 0), the infrared stable fixed point at ii = u* > 
describes quantum liquids of either bosons or ferraions with repulsive interactions which are 
generically universal in the low density limit. Ind = I this fixed point is described by the spinless 
free Fermi gas ('Tonks' gas), for all statistics and spin of the constituent particles, (b) For d > 2 
(e < 0) the infrared unstable fixed point at m = «* < describes the Feshbach resonance which 
obtains for the case of attractive interactions. The relevant perturbation (if — if*) corresponds to the 
the detuning from the resonant interaction. 



Notice the unstable fixed point present for d > 2 and u < 0. Thus accessing 
the fixed point requires fine-tuning of the microscopic couplings. As discussed in 
Refs. lfT2l[m . this fixed point describes a Fermi gas at a Feshbach resonance, where 
the interaction between the fermions is universal. For u <u*, the flow is to m — >^ — oo: 
this corresponds to a strong attractive interaction between the fermions, which then 
bind into tightly bound pairs of bosons, which then Bose condense; this corresponds 
to the so-called 'BEC regime. On the other hand, for u > u*, the flow is to m 
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0, and the weakly interacting fermions then form the Bardeen-Cooper-Schrieffer 
(BCS) superconducting state. 

Note that the fixed point at u = u* for Zf ^ has two relevant directions for d >2. As 
in the other problems considered earlier, one corresponds to the chemical potential 
jj.. The other corresponds to the deviation from the critical point u — u*, and this 
(from ( |36] |) has RG eigenvalue — e=t/ — 2>0. This perturbation corresponds to 
the "detuning" from the Feshbach resonance, v (not to be confused with the symbol 
for the correlation length exponent); we have v °^ u~ u* . Thus we have 

dim[^]=2 , dim[v]=c/-2. (51) 

These two relevant perturbations will have important consequences for the phase 
diagram, as we will see shortly. 

For now, let us understand the physics of the Feshbach resonance better For 
this, it is useful to compute the two body T matrix exactly by summing the graphs 
in Fig. |2] along with a direct interaction first order in uq. The second order term 
was already evaluated for the bosonic case in (l34t for zero external momentum and 
frequency, and has an identical value for the present fermionic case. Here, however, 
we want the off-shell T-matrix, for the case in which the incoming particles have 
momenta ki 2, and frequencies (0\,2- Actual for the simple momentum-independent 
interaction mq, the T matrix depends only upon the sums k = ki+k2 and (o = (Oi + 
a>2, and is independent of the final state of the particles, and the diagrams in Fig.|2] 
form a geometric series. In this manner we obtain 

1 1 



T{k,i(o) Mo 

dQ f d''p 



2n J {2nY [-i[Q. + m) + {p + kf / {2m)) {iQ + p^ / {2m)) 



\ ^ d'^p m ^ r{l-d/2)^j/2 



MO Jo {2%)'^ p^ {Any'l^ 



k^ 

-l(0 + — 
4m 



d/2-l 

(52) 



In c/ = 3, the 5-wave scattering amplitude of the two particles, /q, is related to 
the T-matrix at zero center of mass momentum and frequency k^/m by fo{k) = 
—ml {Q,k^/m)/{47t),and so we obtain 

— L/a — ik 

where the scattering length, a, is given by 

1 47t d^p 4n 

• ' (54) 



a muQ Jo {27t)^ p 



For Mo < 0, we see from ( l54l i that there is a critical value of mo where the scattering 
length diverges and changes sign: this is the Feshbach resonance. We identify this 
critical value with the fixed point m = u* of the RG flow ( |36] |. It is conventional to 
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identify the deviation from the Feshbach resonance by the detuning v 

v = --. (55) 
a 

Note that v °^ u — u*, as claimed earlier. For v > 0, we have weak attractive inter- 
actions, and the scattering length is negative. For v < 0, we have strong attractive 
interactions, and a positive scattering length. Importantly, for v < 0, there is a two- 
particle bound state, whose energy can be deduced from the pole of the scattering 
amplitude; recalling that the reduced mass in the center of mass frame is m/2, we 
obtain the bound state energy, Ei, 

Eb = . (56) 

m 

We can now draw the zero temperature phase diagram liiTI of Zps as a function 
of ji and V, and the result is shown in Fig.|5] 



Superfluid 



BCS 




Fig. 5 Universal phase diagram at zero temperature for the spinful Fenni gas in <:/ = 3 as a function 
of the chemical potential and the detuning v. The vacuum state (shown hatched) has no particles. 
The position of the v < phase boundary is determined by the energy of the two-fermion bound 
state in |56): = — v^l{2m). The density of particles vanishes continuously at the second order 
quantum phase transition boundary of the superfluid phase, which is indicated by the thin continu- 
ous line. The quantum multicritical point at = v = (denoted by the filled circle) controls all the 
universal physics of the dilute spinful Fermi gas near a Feshbach resonance. The universal proper- 
ties of the critical line /x = 0, v > map onto the theory of Section|2] while those of the critical line 
yti = — V^/ (2m), V < map onto the theory of Section|3] This implies that the T >Q crossovers in 
Fig.[T|apply for v > (the "Fenni liquid" region of Fig.[T]now has BCS superconductivity at an 
exponentially small T), while those of Fig.|3]apply for v < 0. 



For V > 0, there is no bound state, and so no fermions are present for /i < 0. 
Ai n —Q, we have an onset of non-zero fermion density, just as in the other sec- 
tions. These fermions experience a weak attractive interaction, and so experience 
the Cooper instability once there is a finite density of fermions for > 0. So the 
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ground state for > is a paired Bardeen-Cooper-Schrieffer (BCS) superfluid, as 
indicated in Fig. |5] For small negative scattering lengths, the BCS state modifies 
the fermion state only near the Fermi level. Consequently as /i \ (specifically 
for < V^/m), we can neglect the pairing in computing the fermion density. We 
therefore conclude that the universal critical properties of the line jj. ~Q, v > map 
precisely on to two copies (for the spin degeneracy) of the non-interacting fermion 
model Zf studied in Section|2] In particular the T > properties for v > will map 
onto the crossovers in Fig. [T] The only change is that the BCS pairing instability 
will appear below an exponentially small T in the "Fermi liquid" regime. However, 
the scaling functions for the density as a function of jx/T will remain unchanged. 

For V < 0, the situation changes dramatically. Because of the presence of the 
bound state ( |56] ), it will pay to introduce fermions even for ji < 0. The chemical 
potential for a fermion pair is 2/i, and so the threshold for having a non-zero density 
of paired fermions is jj. = Ef,/2. This leads to the phase boundary shown in Fig. |5] 
at = — V^/ (2m). Just above the phase boundary, the density of fermion pairs in 
small, and so these can be treated as canonical bosons. Computations of the inter- 
actions between these bosons [11] show that they are repulsive. Therefore we map 
their dynamics onto those of the dilute Bose gas studied in Section[3] Thus the uni- 
versal properties of the critical line pL ~ — v^/(2m) are equivalent to those of Zg. 
Specifically, this means that the T > properties across this critical line map onto 
those of Fig. [3] 

Thus we reach the interesting conclusion that the Feshbach resonance at /i = v = 
is a multicritical point separating the density onset transitions of Zp (Section |2]i 
and Zb (Section [5]). This conclusion can be used to sketch the T > extension of 
Fig-El on either side of the v = line. 

We now need a practical method of computing universal properties of Zps near 
the /i = V = fixed point, including its crossovers into the regimes described by Zf 
and Zb- The fixed point ( [36l ) of Zfs provides an expansion of the critical theory in 
the powers of e = 2 — d. However, observe from Fig. |4] the flow for m < m* is to 
u The latter flow describes the crossover into the dilute Bose gas theory, Zb, 

and so this cannot be controlled by the 2 — d expansion. The following subsections 
will propose two alternative analyses of the Feshbach resonant fixed point which 
will address this difficulty. 



4.1 The Fermi-Bose Model 

One successful approach is to promote the two fermion bound state in ( [56] l to a 
canonical boson field Wb- This boson should also be able to mix with the scattering 
states of two fermions. We are therefore led to consider the following model 

Zfb = / WFt(jC,T)W/r;(A;,T)WB(jt:,T)exp ( - / dTd^'x^FIi] , 
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+ 5 + ^ I VfFal' + (5 - 2li)\^B? 

- Ao (fB*ff t-f^F; + 'i'sfnfn) • (57) 

Here we have taken the bosons to have mass 2m, because that is the expected mass 
of the two-fermion bound state by Galilean invariance. We have omitted numerous 
possible quartic terms between the bosons and fermions above, and these will turn 
out to be irrelevant in the analysis below. 
The conserved U(l) charge for Zfg is 

Q = w;^Wf^ + f;^ff ; + 2%Wb, (58) 

and so Zfb is in the class of models being studied here. The factor of 2 in 
accounts for the 2/i chemical potential for the bosons in ( |57| |. For /i sufficiently 
negative it is clear that ZfB will have neither fermions nor bosons present, and so 
(2) = 0- Conversely for positive /i, we expect {Q) ^ 0, indicating a transition as a 
function of increasing /i. Furthermore, for 5 large and positive, the Q density will 
be primarily fermions, while for 5 negative the Q density will be mainly bosons; 
thus we expect a Feshbach resonance at intermediate values of 5, which then plays 
the role of detuning parameter 

We have thus argued that the phase diagram of ZfB as a function of pL and 5 is 
qualitatively similar to that in Fig. |5] with a Feshbach resonant multicritical point 
near the center. The main claim of this section is that the universal properties of Zfb 
and Zfs are identical near this multicritical point lfTn[T2l . Thus, in a strong sense, 
the theories Zfb and Zps are equivalent. Unlike the equivalence between Zg and Zf, 
which held only 'md= 1, the present equivalence applies for d>2. 

We will establish the equivalence by an exact RG analysis of the zero density 
critical theory. We scale the spacetime co-ordinates and the fermion field as in (|6]l, 
but allow an anomalous dimension r\h for the boson field relative to (|32): 

/ -t 
X = xe , 

t' = xe-^\ 

= Aoe(4-^-')*)^V2 (59) 

where, as before, we have z = 2. At tree level, the theory Zfb with /i = 5 = is 
invariant under the transformations in (|59] with 7]^ — 0. At this level, we see that 
the coupling Aq is relevant for d <4, and so we will have to consider the influence 
of Aq. This also suggests that we may be able to obtain a controlled expansion in 
powers of {4 — d). 
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Upon considering corrections in powers of Aq in the critical theory, it is not dif- 
ficuh to show that there is a non-trivial contribution from only a single Feynman 
diagram: this is the self -energy diagram for which is shown in Fig.|6] All other 



Fig. 6 Feynman diagram contributing to the RG. The dark triangle is the Ao vertex, the full line is 
the propagator, and the dashed line is the propagator 



diagrams vanish in the zero density theory, for reasons similar to those discussed 
for Zb below ( |38] |. This diagram is closely related to the integrals in the T-matrix 
computation in ( |52] |. and leads to the following contribution to the boson self energy 

EB{k,ico) 

dQ d''p 1 1 



2n JAe-i- {2kY' {-i{Q. + (0) + {p + kf/{2m)) {iQ+p^/{2m)) 
^ d''p 1 
'Ae-t {2nY {-ico+{p + kY/{2m)+p^/{2m)) 
^ d''p m f 4\\ /-^ d'^p 

The first term is a constant that can absorbed into a redefinition of 5. For the first 
time, we see above a special role for the spatial dimension d — 4, where the momen- 
tum integral is logarithmic. Our computations below will turn to be an expansion 
in powers of (4 — d), and so we will evaluate the numerical prefactors in (l60t with 
d = 4. The result turns out to be correct to all orders in (4 — c/), but to see this explic- 
itly we need to use a proper Galilean-invariant cutoff in a field theoretic approach 
ifTTl . The simple momentum shell method being used here preserves Galilean in- 
variance only in d — 4. 

With the above reasoning, we see that the second term in the boson self-energy in 
(|60] | can be absorbed into a rescaling of the boson field under the RG. We therefore 
find a non-zero anomalous dimension 

r7fo = A^ (61) 

where we have absorbed phase space factors into the coupling A by 

Ao = =A. (62) 
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With this anomalous dimension, we use ( |59] l to obtain the exact RG equation for 

A: 

dX {A-d) A3 

dI = ^^"T- ^^^^ 



For < 4, this flow has a stable fixed point at A = A* = \J{A — d). The central claim 
of this subsection is that the theory Zps at this fixed point is identical to the theory 
Zfs at the fixed point u — u* for 2 <d <A. 

Before we establish this claim, note that at the fixed point, we obtain the exact 
result for the anomalous dimension of the the boson field 



rih^4-d. (64) 

Let us now consider the spectrum of relevant perturbations to the A = A* fixed 
point. As befits a Feshbach resonant fixed point, there are 2 relevant perturbations 
in Zfb, the detuning parameter 8 and the chemical potential fi. Apart from the tree 
level rescalings, at one loop we have the diagram shown in Fig.|7] This diagram has 



Fig. 7 Feymnan diagram for the mixing between the renormalization of the W^ff and "f^ 
operators. The filled circle is the "f^ source. Other notation is as in Fig.|6] 

a ^p^^fa source, and it renormalizes the co-efficient of <P^<P; it evaluates to 
dQ /-^ d''p 1 



2 d'p 



Combining ( |65] | with the tree-level rescalings, we obtain the RG flow equations 

^(5-2m) = (2-77i)(5-2Ai)-2A2^, (66) 
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where the last term arises from (|65t . With the value of rjt, in (|6TT i. the second equa- 
tion simplifies to 

f = (2-A^)«. (67, 

Thus we see that ji and 8 are actually eigen-perturbations of the fixed point at A = 
A*, and their scaling dimensions are 

dim[ju]=2 , dim[5] ==t/-2. (68) 

Note that these eigenvalues coincide with those of Zps in ( BTI ). with 5 identified as 
proportional to the detuning v. This, along with the symmetries of Q conservation 
and Galilean invariance, establishes the equivalence of the fixed points of ZfB and 

The utility of the present ZfB formulation is that it can provide a description of 
universal properties of the unitary Fermi gas in (i = 3 via an expansion in {4~ d). 
Further details of explicit computations can be found in Ref. lfT2l . 



4.2 Large N expansion 

We now return to the model Zps in ©, and examine it in the limit of a large number 
of spin components ifTTllZTl . We also use the structure of the large perturbation 
theory to obtain exact results relating different experimental observable of the uni- 
tary Fermi gas. 

The basic idea of the large expansion is to endow the fermion with an addi- 
tional flavor index a = 1 . . .N/2, to the fermion field is ^pan, where we continue to 
have (7 =t, i- Then, we write Zfs as 

Zf, = /WF„„(x,T)exp(-/J/^^/T/flf''x^f,) , 

■^F^ = I VfFa« I' " ^aal' (69) 

where there is implied sum over a,b^l .. .N/2. The case of interest has N ~2, but 
we will consider the limit of large even A^, where the problem becomes tractable. 

As written, there is an evident 0{N/2) symmetry in Zfj corresponding to rota- 
tions in flavor space. In addition, there is U(l) symmetry associated with Q conser- 
vation, and a SU(2) spin rotation symmetry. Actually, the spin and flavor symmetry 
combine to make the global symmetry U( 1 ) x Sp(A^), but we will not make much use 
of this interesting observation. 

The large expansion proceeds by decoupling the quartic term in (|69^ by a 
Hubbard-Stratanovich transformation. For this we introduce a complex bosonic field 
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Wb(x,t) and write 

Zf, = / Wfa„(x,T)WB(x,T)exp J^'' dT jd^X^Fs) , 

= 'f'fV, + 2k |Vff aal' - Miff a„P (70) 



/|Mo| 

Here, and below, we assume hq < 0, which is necessary for being near the Feshbach 
resonance. Note that couples to the fermions just like the boson field in the Bose- 
Fernii model in (|57] |. which is the reason for choosing this notation. If we perform 
the integral over ffe in ( fTOl i. we recover (|69), as required. For the large expansion, 
we have to integrate over f^aa first and obtain an effective action for Because 
the action in (iTOl l is Gaussian in the ^Faa, the integration over the fermion field 
involves evaluation of a functional determinant, and has the schematic form 

J>. = y"wB(x,T)exp(-A^^eff[fs(x,T)]), (71) 

where ^eff is the logarithm of the fermion determinant of a single flavor. The key 
point is that the only dependence is in the prefactor in (iTll . and so the theory of 
can controlled in powers of 

We can expand ,5^^if in powers of ffe: the p'th term has a fermion loop with p 
external "ffe insertions. Details can be found in Refs. ifTTl l27l . Here, we only note 
that the expansion to quadratic order at /i = 5 = T = 0, in which case the co-efficient 
is precisely the inverse of the fermion T-matrix in ( |52] |: 

- 4/g^^m(*,a,)P + ,„ (72, 

Given S'eff, we then have to find its saddle point with respect to f^. At T = 0, we 
will find the optimal saddle point at a ffe 7^ in the region of Fig. |5]with a non-zero 
density: this means that the ground state is always a superfluid of fermion pairs. The 
traditional expansion about this saddle point yields the expansion, and many 
experimental observables have been computed in this manner lfmi27ll28l . 

We conclude our discussion of the unitary Fermi gas by deriving an exact rela- 
tionship between the total energy, E, and the momentum distribution function, n{k), 
of the fermions Il25ll26l . We will do this using the structure of the large expan- 
sion. However, we will drop the flavor index a below, and quote results directly for 
the physical case ofN = 2. As usual, we define the momentum distribution function 
by 

n{k) = {W^Jk,t)^Fa{k,t)), (73) 

with no implied sum over the spin label a. The Hamiltonian of the system in ( |69l ) is 
the sum of kinetic and interaction energies: the kinetic energy is clearly an integral 
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over n{k) and so we can write 



E =2V 



f d''k k^ , , dXnZps 

where V is the system volume, and all the fields are at the same x and f . Now let 
us evaluate the mq derivative using the expression for Zps in ( iTOb ; this leads to 

1=2/ + (x,0'f^^(x,f)). (75) 

y J \1%Y 2m Mo 

Now using the expression (|54] | relating mq to the scattering length a in c/ = 3, we can 
write this expression as 

This is the needed universal expression for the energy, expressed in terms of n{k) 
and the scattering length, and independent of the short distance structure of the 
interactions. 

At this point, it is useful to introduce "Tan's constant" C, defined by 12511261 

C^\imk*n{k). (77) 

The requirement that the momentum integral in (|76) is convergent in the ultraviolet 
implies that the limit in ( ITTI i exists, and further specifies its value 

C = m2(fB*fB). (78) 

We now note that the relationship n{k) {W^Wb) jk^ at large k is also as 

expected from a scaling perspective. We saw in Section l4n that the fermion field 
does not acquire any anomalous dimensions, and has scaling dimension t//2. 
Consequently n(ki) has scaling dimension zero. Next, note that the operator "fg "ffe is 
conjugate to the detuning from the Feshbach critical point; from (|68] | the detuning 
has scaling dimension c/ — 2, and so "fg "ffe has scaling dimension d+z— (c/ — 2) = 4. 
Combining these scaling dimensions, we explain the k^'^ dependence of n(]i). 

It now remains to establish the claimed exact relationship in dTST l as a general 
property of a spinful Fermi gas near unitarity. As a start, we can examine the large 
k limit of n{k) in the BCS mean field theory of the superfluid phase: the reader 
can easily verify that the text-book BCS expressions for n(k) do indeed satisfy (|78] ). 
However, the claim of Refs. ||Tl|24l is that ( iTSl l is exact beyond mean field theory, and 
also holds in the non-superfluid states at non-zero temperatures. A general proof was 
given in Refs. 124], and relied on the operator product expansion (OPE) applied to 
the field theory (iTOl i. The OPE is a general method for describing the short distance 
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and time (or large momentum and frequency) behavior of field theories. Typically, 
in the Feynman graph expansion of a correlator, the large momentum behavior is 
dominated by terms in which the external momenta flow in only a few propagators, 
and the internal momentum integrals can be evaluated after factoring out these fa- 
vored propagators. For the present situation, let us consider the 1 /N correction to 
the fermion Green's function given by the diagram in Fig. |8] Representing the bare 




Fig. 8 Order 1 /N correction to the fermion Green's function. Notation is as in Fig.|6] 



fermion and boson Green's functions by Gf and Gb respectively, Fig[8]evaluates to 

^^—GB{p,a)GF{-k + p,-(0 + a). (79) 

Here Gb is the propagator of the boson action o5^eff specified by ( l72l ). In the limit 
of large k and ft), the internal p and Q. integrals are dominated by p and £1 much 
smaller than k and ft); so we can approximate ( |79] l by 

gI{K(0)Gf{-K-(0) J ^^GBip,£2) 
= Gl{k,(o)GF{-k-a){%^'B). (80) 

This analysis can now be extended to all orders in Among these higher order 
contributions are terms which contribute self energy corrections to the boson propa- 
gator Gb in dSOl i: it is clear that these can be summed to replace the bare Gb in dSOl l 
by the exact Gb- Then the value of (I'ffep) in (ISOl i also becomes the exact value. All 
remaining contributions can be shown fl^ to fall off faster at large k and ft) than 
the terms in dSOl l. So dSOl l is the exact leading contribution to the fermion Green's 
function in the limit of large k and ft) after replacing (jffep) by its exact value. We 
can now integrate (ISOl l over ft) to obtain n{k) at large k. Actually the ft) integral is 
precisely that in (l65T l. which immediately yields the needed relation ( fTSl l. 

Similar analyses can be appUed to determine the the spectral functions of other 
observables ll28l l8l l20l |2T1 HI El l24l . 

Determining of the specific value of Tan's constant requires numerical compu- 
tations in the 1 /N expansion of (TtTI i. From the scaling properties of the Feshbach 
resonant fixed point in = 3, we can deduce the result obeys a scaling form similar 
to (Eli; 
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where <Pc is a dimensionless universal function of its dimensionless arguments; note 
that the arguments represent the axes of Fig. |5] The methods of Refs ifTTllZTl can 
now be appHed to ( |78] l to obtain numerical results for <t>c in the 1 /N expansion. We 
illustrate this method here by determining C to leading order in the l/N expansion 
at /i = V = 0. For this, we need to generalize the action (iTZt for to T > and 
general A^. Using ( |52] i we can modify (|72]i to 

yefi^NTV [ P^[Do{k,(0„)+Diik,(On)]\1'B{k,0)n)\\ (82) 

where Dq is the T = contribution, and Di is the correction at 7" > 0: 



Do{k, (0„) = -—\ -i(0„ + — (83) 
Ion V 4m 

\ rd^p 1 1 

D\{ .®")-2y 8;r3 (gPVlZmr) + j) {-i(o + p'^/{2m) + {p + kf/{2m)y 

We now have to evaluate (f^f^) using the Gaussian action in ( [82] i. It is useful to 
do this by separating the Dq contribution, which allows us to properly deal with the 
large frequency behavior So we can write 



1 1 



DQ,{k,cOn)+Dx{k,co„) Do{k,co„) 



+ DoQ. (84) 



In evaluating Dqo we have to use the usual time-splitting method to ensure that the 
bosons are normal-ordered, and evaluate the frequency summation by analytically 
continuing to the real axis: 

I f d^k ^ e''""'' 
Dqq = — — 5- lim T > 



Nj 87r3,,^o ^Do(fc,co„) 
167r f d^k r da 1 



A?ot3/2 J 87r3 n {e^lT - 1) ^ a - / {Am) 

= (85) 

The frequency summation in (|84] | can be evaluated directly on the imaginary fre- 
quency axis: the series is convergent at large o,,, and is easily evaluated by a direct 
numerical summation. Numerical evaulation of (l84t now yields 

C={2mT)^{^^ + ff{l/N')) (86) 



at /i = V = 0. 
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